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SUBMERGED | BROAD. ) CRESTED WEIR 


Richard A. . Smith, 1M. 


| 


2 treats the calibration, or correlation of variables, even though the weir 7 
rae submerged and tailwater phenomena control the discharge.  ‘Pre- 

_ vious writings treat the calibration for the condition of free over-fall, but for 

{ _ the submerged condition, the discharge coefficients were shown to vary over _ 
a wide range in unpredictable fashion. This paper shows that for fr free over- 

fall the discharge is a function of the head on the weir, but for the . submerged > 
"condition the discharge is a function of the degree of submergence. 
_ calibration, the computed discharge, when compared with the observed dis- 
3 “charge, § showed a probable deviation of approximately 2 percent. ine» 


Studies of floods in Louisiana entail the connie: of discharge over ~~ 
inundated levees or roadbeds. Available formulas appeared inadequate a 
_ computing the discharge over the levee or roadbed when considered as a 
broad crested weir. _ With available weir formulas, their accuracy was vitally: 
_ dependent on empirical coefficients that were little understood in many in- . 
stances. | In application, a value was somewhat arbitrarily assumed without = 


fi 
= knowing whether the chosen value wa was even . approximately correct. 


At the request of the State ai Louisiana, Department of Public Works, hy- 
- draulic model studies were instituted. The purpose was to calibrate the ey: 
1 inundated roadbed, acting as a weir, for two conditions of flow: first, when 
x. Note: Discussion open until August 1, 1959. To extend the closing date one month, a 

_ written request must be filed with the Executive Secretary, ASCE. Paper 1960 is 

= of the copyrighted Journal of the Hydraulics Division, Proceedings of the a 


American Society of Civil Engineers, Vol. 85, No. HY 3, March, 1959. 
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controlled the. discharge; when the tail-v ~water | was igh, 


‘yield elevations of the v weir - crest, head- water, ‘and tail-water surface a 


In the analysis of data, the epecific objectives were threefold 


a) Determine a criterion by which it could be ascertained when control of 
regimen passed from the submerged weir to tail- 


b) Calibrate. the head- discharge relationship for the condition when c con- 


trol of the regimen was atthe weir, 
Correlate t the variables, head-water elevation, tail-water elevation, 
depth of flow over the weir crest, and rate of discharge, for the condi- . 
when control of the regimen was downstream from the 


slope of the upstream and faces. . Therefore, model studies were 
conducted for two models having different slopes. One had slopes’ of 1 on 1; =e 
_ the other had slopes of 1 vertical on 2 horizontal, with the upstream slope pei ' 
" being the same as the downstream slope it in either ¢ case. Other itera sey: 
- mon to both models, were as follows: Length of weir r normal t to the direction | 
of am, 4. 95 feet; width of weir crest parallel to the direction of flow, 0.5 a 
A 
roughness was established by sprinkling the models v with masons sand shortly 
after they had been painted. A diagram of the test apparatus is shown ge ‘ 


Head- ‘water and tail-water elevations were controlled by a supply valve 
oh and tail gate, respectively. _Head-water and tail- water elevations were ob- 

served by means of piezometric connections, 1/2 inch diameter, to stilling 
wells equipped with hook gauges. The depth of flow over the weir crest was | 
_ observed by means s of a point gauge located slightly downstream from ST 
- centerline. This point gauge was read when the fluctuating water surface ‘ 
was | above the point and below the point for approximately equal intervals of 


et = (Piezometric connections were e not deemed feasible at the weir crest — 


observed means of a weighing tank and electric ‘stop- clock, “which read 
“directly to one-tenth of a second. | ‘The model was operated throughout the %. 

range of discharge capacity available, , up to 1. 0 cfs on the weir five feet long, © 3 
or 0.2 cfs per foot of weir length. At each rate of discharge, the tailwater Ls e 
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epth of flow. over the weir crest | 
eleva i mn above the weir crest 


4 
“= | 
i 
| 
q 
re 
J 
| 
| 
Symbols 
es- 


‘Tables 1 and 2 list data observed for the sinitate with slopes. oft the weir 


faces 1 on 1 and 1 on 2 


; ‘ing the circumstances when control of the regimen passed from the 1 sub- f - 


merged weir to tail-water conditions, or conversely: _ 


} mi When the weir was controlling the regimen, the tail- water elevation 


b) When control of the regimen was in the tail-water conditions down- 


a stream from the weir, each change in tail-water elevation produceda 
corresponding change in the depth of flow over the weir crest and in 


head-water elevation, with the discharge remaining constant. 


c) Between successive runs, when control passed from the weir to tail- 
ha a water conditions, or conversely, this change in regimen was reflected 


by marked changes in the Froude number of flow over the weir crest. We 


to the tailwater, and conversely, occurred in the region of Froude num- 

- bers approximating 0. -78 for flow over the weir crest. | The higher num- 
bers correspond to the control at the 1 weir, while numbers less than — " 


0.78 correspond to control in the tail- water conditions. 
vie 4, 


Coefficients, Weir Controlling 
From Figure 2, the head- -discharge ‘relationships f for the models were com- 


variable, ‘this | yielded the following 


the model lonl slopes and, 
64H 
having of 1 on 2. A graph of 
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CRESTED WEIRS- 


OBSERVED DATA FROM ‘DISCHARGE’ OVER MODEL OF OF ROADBED 
-HAVING SIDE SLOPES OF 1ON1 


CREST TAILMATER DISCHARGE OBSERVED 
«DEPTH ‘PSR FOOT COSPPICIENT NUMBER 


0. 034 025 = 000 0183 
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“TABLE 1 (CONTINUED) 


‘DEPTH BLEVATION FOOT COBFF ICLENT 
0. 096 084 0. 077 a 0839 2. 
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‘DEPTH BLEVATION FOOT NUMBER 


0.107 0.085 0.089 1417) 
6. 054 1405 
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(©) Observed besttictont is markedly different from computed mean, but 
the observation is in the range of ‘Small heads where large 


€rrors may be anticipated. 


The weir is drowned out, and the head-discharge relationship ts 


controlled by tail-water phenomena. 


Coefficients, Tail -water 
_ In Figure 3 is shown graphs of the degree of pubmergence and relative e 


depth at the weir crest as related to the Froude number. The equations of 
‘the curves and the graphs were derived by the method of least squares. In 
‘ order to avoid undue confusion of the graph, the individual data points are ; 


omitted. 
of 


For the model on the weir faces of 1 on the 


submergence is related to the Froude me equation, 
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TAILAATER- DISCHARGE OBSBRVED 
DEPTH PER ‘FOOT  COBPPICIENT 
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TABLE 2 


HEAD DEPTH ELEVATION PER FOOT CORFFICIENT NUMBER 
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) Observed coefficient is markedly different mean, but” 

the observation is in the Senge of Small heads where 


_ The relative dines at one weir crest is related to the Froude nun number by 
‘the equation, obtained by least squares analysis: 


— 
RUN HEAD 
| 
ol ae 
- 
the resulting equation is as follows: 
p <2 By comparing this equation with the 
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HEAD -DISCHARGE 


ACTING AS WEIRS. 
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the on the weir faces of 1 on 2, the of sub- 
is to Froude number by the equation, obtained by least 


a 
Rare 
£455 
4 
a | = 
: q are 
4 
4 | 
4 


SUBMERGENCE 


OF 


EE 


& 


a 


(RELATIONSHIP OF THE 
FROUDE N UMBER TO DEGREE _ 


TH, FOR MODEL. S OF ROA 
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equations (5) and ( ) are with the continuity equatio! 
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g equation is as follows: 


this equation with the weir 


2. 66 - 


Thus it is seen in equations (4): and (8) ‘that the coefficients of discharge 
: “a either model is a function of the degree of submergence (z/H). It should ‘~ 
be noted that this correlation is for the condition of the weir drowned out < a 
- and control of the regimen is in the tail-water. A graph of these coefficients i. 
is shown in Figure 4, representing equations (4) and (8), above. = 


In compe 


aring observed data with the it was noted 
that large deviations occurred in 2 runs out of a total of 53 for one model and 
4 in 4 runs out of a total of 46 for the other model. In each instance of large | LA 
. deviation, the total head from head- -water to tail -water was 0.018 feet or less. is 
- Figure 5 shows that for smaller heads, an error in observed head will pro- aa 


duce a larger percent error in computed discharge. _ Therefore, in arriving 
> a at the probable deviation, the . larger deviations v were omitted from the a <4 
analysis. Those runs omitted are indicated in Tables 1 and 2, sca parentheses: > 
enclosing th the numerical value of in Figure 4. 
Errors i in Discharge Relative to in Head ne 
The following argument will show the relative deviation in| 
caused a corresponding deviation in observed head. weir 
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4 DEGREE OFS: UBMERGENCE 
TAIL- WATER CONTROLLING — 


OF ROADBEDS 
WITH CONTRO L EITHER A 
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For small, finite increments, this” may expressed as follows: 


“The percent | deviation in computed discharge is 1. 5 times the percent Bice 94 
deviation in observed head. This relationship is shown ail 
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of Observed Data Fr From Computed 


The! laws or probability were applied to the. observed 


- with mean values computed | by the weir formula and using coefficients from 
Figure 4. comparison n indicated that the observed data agreed with the 


computed mean vanes with a probable deviation of a plus or minus 1.4 
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arise, in n open channel flow may be solved without tr 


@ 


— 
~ Most engineers will cheerfully admit their inability to solve equations be- 
~ yond quadratics by any method other than “cut-and- -try.” Special methods for 
a cubic equations have been devised for particular equations of frequent R: 
_ occurrence, and such a method is presented here for an equation that arises — 


“Although this diagram is reproduced in practically every fluid mechanics _ 
: textbook on the market, (1) none of them explain the use of the diagram in 
practical calculations. Figure 1 shows a specific energy diagram plotted to 
F . scale in dimensionless form. Actually, a diagram cannot be read with suffi- a 
cient precision for many purposes and so a giving me same data is also 
a rectangular specific energy at any eection 


= “Discussion open extend the closing one month, 
written request must be filed with the Executive Secretary, ASCE. Paper 1961 is 
_ part of the copyrighted Journal of the Hydraulics Division, —— of the oe 
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symbols have the meanings: 
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EQUATIONS 
Values of for values sof =! 4.99 are given in Table 


1 was plotted from Table 1 


fl by appropriate manipulation reducing it to the form of Eq. 3. Whenthis — 
ae _is done, « one side or the other | of Eq. 3 can be evaluated from given data, and — 


1. Finally, multiplication of tabular ‘a 
values by de will yield a numerical 
_ If it is necessary to include an energy loss between two cross sections a 
the table or graph may be used by making a slight modification of the theory | 


above. Say that between section and the 


section there is a loss equal to K 3° . (The subscript 2 refers to the down- 
t 


stream section. ) Bernoulli’s s theorem written aiid the upstream | se 
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m* Te and 


Division n of Eq. 5 by Bq. gives 
al 


This equation is s of the same form as Eq. 3 and when either side is nice ’ 
the other side may | found from Table 1 or 
Llustrative Probl 4 
“The problems would lead to cubic equation if accompanying 


assure accuracy of the decimal 
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g/d, F 


0.02 6 0. 0. 0. 
5000 1250. 555. 312. 5 200.1 138.9 102.1 78.21 61.82 
43 34, 84 29.71 25.65 22.37 19, 69 17.47 ‘As. 61 
.2 12.70 11.55 10.55 9.682 8. 920 8.250 7.656 7. 129 6.658 
856 5. 5. 203 4.921 4.665 4.432 218 022 3. 3 3. 
.525 3.384 3.254 3.134 3. 023 2.919 2.823 2.734 «2.650 2. 
.500 2. 432 2.369 2.310 2. 255 2.203 2.154 2.109 2.066 2. 
.989 1.954 i. 921 1,890 1.861 1. 833 1. 784 A 761 1. 
1.702 1.684 1.653 1.639 1.626 
"572 1. 564 1.556 1.549 1. 542 1.536 1. 526 
500 .502 1. 504 507 1.509 1.511 
516 1. 519 1.522 1.525 1.528 1.532 1.535 1.539 
-552 1.556 1.560 +565 1.570 1.575 1.580 21,585 1.590 
.601 1.607 1.613 1.618 1. 624 1.630 1.636 1.642 i 
.662 1.668 -681 1.688 1. 695 1. ..708 
(729 1.736 1.744 1.751 1. 1. 773 1.780 
.803 1.810 1.818 1.834 1. 841 "349 1 1.865 
881 1.889 1.897 -905 1.913 1. 921 1 11930 1. 938 1.946 
"963 1. 971 1. 979 1.988 1.996 2.004 2.013 2.022 2.030 
.047 2.056 2.064 2.073 2.082 2.090 2.099 2.108 2.116. ; 
2, 142 2. 2. 160 2.169 2.178 2.187 2. 196 2. 204 
222 2.231 2.908 2. 249 2.258 2.267 285 2.294 
2. 322 2. 330 349 2.358 2.367 2 376 2. 385 
2. 431 2.440 2.450 468 2. 
-496 2.505 2.515 2.524 2.533 2.543 2.552 2.561 > 
589 2. 599 2.608 2.618 2.627 2.636 2. 646 2.655 2. 
2, 693 2.20 2.7 2 2. 750 2. 
.778 2.788 2. 807 2.816 2,826 2.835 2.845 2. 
873 2. 883 (2. 2.902 (2,912 2. 921 2. 931 2. 940 
998 3.008 3. 3.027 2.036 
065 3 3.075 3. 3.094 3. 3.113 3.123 3.133 
162 3.171 3.181 3.191 3.200 3.210 3.220 3.229 3. 
3.268 3. 3.288 3.297, 31307, 3.317 3.326 3. 
365 3. 3.385 3. 395, 3.404 3.414 3.424 3. 
453 3.463 3.472 3.482 3. 492 3.502 3.512 3.521 3. 
(551. 3.560 3. .580 3.590 3. 600 3.609 3.619 3. 
.648 3.658 3. 678 3.688 3.697 3.707 3. 717 3. 
.746 3. 756 3.766 3.776 3, 786 3.795. 3. -805 3.815 3. 
3.844 3.854 3. .874 3.884 3.894 3. 903 3.913 3. 
943 3.953 3 .972 982 3.992 002 4.012 4. 
031 4.041 4 Osi 4.061 4.071 4.080 4.090. 100 4 110 4. 
.130 4, 140 4, 150 4. 159 4.169 4.179 4.189 4. 199 4. 209 4 
.228 4.238 4.248 4.258 4. 268 4.278 4.288 4.297 4. 307 4. 
4.337 4.347 4.357 4.366 4.376 4.386 4.396 4.406 4, 
+426 4 436 4. 446 4.456 4.465 4. 475 4. 485 4.495 4.505 4. 
4.525 4.535 4.544 4 4.554 4.564 4.574 4.584 4.594 4.604 4 
4.634 4.653 4.663 4.673 4.683 4.693 4. 
.723 4.732 4. a 4.752 4.762 4.772 4.782 4.792 4.802 4 
4, 832 4, 842 4,851 4.861 4.871 4.881 4.901 4. 
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4 Example 1 


az of 7 feet be carried by a rectangular channel 9 feet eee et - — 
th the given data the following: 
3.945 


3. 9450 


xi 674. Multiply these v values by de to obtain 
6. 198 feet (tranquil quil flow) 


feet wide and deep carries 300 cubic feet p per second. 
‘ A transition reduces the width to 9 feet with no change in ‘bottom elevation. 
_ Find the depth in the narrow section. 
From the given data compute the 


‘Enter 1 with this of — and obtain 

s-1 = 


In the narrower section 


losses are neglected bes Es. $8, and 


q 

F 

and 
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ita a rectangular channel 8 feet wide carries a flew of 250 cubic feet per 
ata depth of 5 feet, change in surface e élevation will be produced 


by a local rise in asks level of 1/2 foot? 


nter Table | 1 with this v value of — to obtain 


‘Divide both sides of this equation by ae - pearnange the terms to obtain 

All quantities on the right of this equation are mews. 


: 
‘Enter Table with this value an and obtain 1. 371 from w which 
the surface will drop 0. 22 ft. 


reservoir. If the entrance loss is 0.2 =— find the water depth in the flum 

a flow of 70 cubic feet per second per foot. | 
Write Beracellt’s sevens between the surface and the channel 


q and combine the last two terms. 


lp = 1.808 x 3.257 = 5.889 feet, — 
aa 
é 
j 
a ag a ye 
| 


m: The e result 


m* 
In 
This is the same form as Equation 3 and Table 1 may be used with dn 


¥ ‘rej lacin do. 


Enter the table with this value ot = 1.558 


558 x 5. 8.83 ft. 
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York, N. Y., problem 161. 
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RELATIONSHIPS BETWEEN PIPE RESISTANCE FORMULAS 
Walter L. Moore,! A.M. ASCE 
The relationship between modern concepts of pipe resistance and the | older 
empirical formulas is clarified, and a simple procedure developed to derive 
an exponential formula applicable to a known range of flow conditions. The a 
analysis makes apparent the range of flow conditions for which any particu-— _ 
ag 
lar exponential formula is applicable. Limitations of the equivalent pipe con-_ a 
cept are discussed. Data on resistance measurements in water mains in 
service indicate that the head loss varies with the discharge to a power near- 
ly equal to 2.0 rather than 1.85 as commonly assumed. 


_ INTRODUCTION 


Methods for predicting the resistance encountered by a fluid flowing 

‘tarough a pipe have long been of concern to mathematicians and engineers. 
; The problem has | been recognized as an important one in the ‘design of water 

supply lines and water distribution systems as well as in process plants of 
_ various types. Recently the use of long pipelines for transporting large _ 
_ volumes of fluids” has renewed ir interest in problems of pipe resistance. _— 

_ Over the years, ‘different approaches have been made to the problem of | 
wf predicting pipe resistance and a multitude of equations | and graphs have been 

_ devised for evaluating the resistance. In recent years, , the literature on the 

over a presents semi- -rational leading to relations 


limitations. Perhaps this should not be surprising because so little has been ns 


written to explain the interrelationship between the various methods of 
Discussion open until August 1, 1959. To extend the date one month, 
__- written request must be filed with ‘the Executive Secretary, ASCE. Paper 1962 is — 
sats ~ part of the copyrighted Journal of the Hydraulics Division, Proceedings of he 
i= American Society Vol. 85, No. HY 3, March, 
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evaluating resistance paper attempts to clarify this interrelationship | 
- ‘showing exactly how the older empirical equations may be regarded as spe- 


cial cases of the more recent general resistance relations. Itis hoped = 


Cent information as it is presented in the general resistance diagram. “This ie = 


Development of Pip e Resistance Relations 
A brief review of some selected major steps in the development of pipe > 
Zz resistance relations will help to clarify the resistance relations to which 
' om is made. In Table I, some selected major steps in the develop- | 
of pipe resistance relations are shown giving the approximate date 
the development, the name of the investigator commonly associated with it, 
the resistance relation, and the conditions for its application. eee 
Z Some of the resistance relations, namely the Chezy, Manning, and Hazen- = 
Williams equations, are customarily in a form convenient for solving for the - 
velocity V. In Table I, these are also written in a form in which the head 
loss hg is expressed in terms of the other variables for the the 


 glope S$ as = In this form, they | are more comparable to the | common form Me 


of the Darcy-Weisbach equation. The relations listed as previous to 1930- 31 “te r 
in Table I are all based on empirical studies with the flow of water. - They “= 


8 


where C isa a coefficient to > account for the pipe material, L the nasal of the 


pipe, R the hydraulic radius, V the mean velocity, m are expo- 


where K isa coefficient which upon the ma- 
, terial of the pipe and Q is the discharge. Equations of this type are sated i 
xponential equations because the head loss | is expressed in terms of the pipe 


The modern understanding of pipe resistance began about 1930 with the — 

a Prandtl-Karman relations which were | supported by the experiments of Asa ad 

-Nikuradse, Colebrook, and White. These modern relations are based on an [ 

| understanding of the variation of the resistance coefficient f in the je fluid + 
4 _ Weisbach equation. . For the steady uniform flow of an incompressible fluid - 
“the complete picture is well summarized in a general resistance diagram ‘e 

such as proposed by Rouse and by Moody in slightly different forms. The © 

values of f from the general diagram are substituted in the Darcy 
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Selected Steps in the of Pipe Resistance Relations 


VRS Flow of water in open channels. 


ite for Laminar Flow in pipes. 


Freeman Exhaustive tests and Flow of water in commercial pipes. 


Prandtl- 4 Turbulent flow of any in 


Colebrook 2 2.52 , e Commercial turbulent flow, 


White "smooth," "rough," ont transition 
*3.D 


 COmmexcial pipe: Laminar and tur-— 
(lige Trims bulent flow, smooth, rough, and 


= orf "Moody Diagram" - for commercial 
¢ (gs ~ pipe: Laminar and turbulent flow, 
R’D smooth, rough, and transition 


= | rough 
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— 3 ____ Each of the two types of formulas, the exponential type and the one based BP 


on the general resistance has | some as well as disad- 
vantages. With the exponential type, the calculations are convenient; and for 
this reason it is widely used by many practicing engineers. The calculations © 
direct, avoiding the need for trial-and-error solutions. The reliability 


the calculated resistance is questionable, however, because of the uncer - 


3 in selecting the applicable equation from the very large number of ex- 


ponential formulas that have been proposed. As late as 1951, Blair(10) refers Ss. 
to the multitude of formulas already proposed and proceeds to suggest four — 
F new exponential formulas designed to cover a wide range of pipe materials. 


fluid they are necessarily limited to conditions where 


The modern general resistance ‘diagram, ‘however, accounts | ina 
manner for all of the factors which influence pipe resistance. A clear under - 
_ Standing o of the mechanism of pipe resistance has come from this approach, 
“4 and the general resistance diagram clearly shows the effect of the various © 
i factors entering into the problem. The improved knowledge incorporated in 
‘ the general resistance diagram, however, has not received as wide usage > 
~ among practicing engineers as it deserves. a ‘This is probably largely due to 
its inconvenience as compared with the simpler exponential formulas. ‘With 
_ the general resistance diagram in its common form, trial-and-error solutions — 
are required for finding the rate of flow or the size ‘of pipe. ie oily ani 
Because of the convenience of the exponential formulas and their wide- 
- spread usage, it is } desirable that the Teeaiamae "ae between them and the 


For a limited range 2 of numbers and a fixed value of 5 

of log f versus log NR may be approximated a line 


as 


log f versus log NR on the If this is substituted into od 
Me oi Darcy- -Weisbach equation and the greeyren number is expressed in “gr of 


a 
a 
| = 
4 
im 
Apain in these 10ormuias ere ala Cieveh new CApOonential 
ae 4 formulas (divided into four groups) were proposed to represent a wide range ri 2 
A 
possible to benefit from the convenience of the exponential formulas as well 4 
as the completeness of the general resistance diagram. The interrelationship 
— _ between exponential formulas and the general resistance diagram has not | Ae ( 
previously been presented in complete form. 
ted incomplete form. 
4 
— 


PIPES 
The ond loss may then be expressed - 


Int this expression both the coefficient Co and th the exponent m will ~ tea al 


4 tions of the epee: number and the relative roughness of the | pipe =. 
when the head loss is expressed in the form of equation (7), the coefficient E> 
Ko depends upon the pipe characteristics of length, diameter, and relative q an. a 
"roughness as well as the Reynolds number and fluid viscosity (see equation (8)). 


_ The yey m may be conveniently determined from the Moody ae, 


: Be of etn m may be drawn on the Moody diagram by connecting points 
of constant slope. This has been done in the modified Moody diagram of Fig. 1. 2 
modified 
las have been proposed. It is ; evident that to completely r represent the r resis- - 
 temee problem an infinite family of exponential formulas is needed with the © 
exponent m varying continuously from at minimum value of about 1 1.7 for 
- smooth pipes at a Reynolds number of about 4000 to a maximum value of 2 — 
for a relatively rough pipe at large Reynolds numbers. For a given value of _ 


use of the proper exponent m is when the head loss is to be 
3 determined for different discharges or when a change of head loss is to be hy 


related to a change in discharge as in the Hardy Cross network analysis. Dif- 
ferentiation of the equation hg = eld — = that the 


Any of the common exponential equations such as the Hazen-Williams, — 
; Manning, or Darcy-Weisbach (with constant f) “may be made to yield the jan 
correct head loss at a given discharge by proper choice of C, n, or f. “How- 
_ ever, the head loss at other discharges as well as the relation between 4 % 
_ change in head loss and a change in discharge will be incorrect unless the yin 
: exponent m of the selected formula corresponds to the combination of rela-— 
_ tive roughness and Reynolds number for the pipe under consideration. The | 
a modified Moody diagram gives an easy check on the proper m value for ‘any 
* a combination of relative roughness and Reynolds number. Thus the Hazen- v 
_ Williams formula with its exponent of 1.85 corresponds to a combination “a ae 7 
relative roughness and Reynolds number lying along the line m= = = 1. 85 on the ‘ran 
Moody diagram. The Manning equation or the Darcy- -Weisbach (with constant a 


9 corresponds to conditions lying on or to the right of the line m = 2. 0. 


"exponential formula for a reasonable range of conditions which will accurately — 
reflect the resistance characteristics as embodied in Moody |The 


steps have been found most convenient: 


= 
4 
> 
____relative change in head loss is_m times the relative change in discharge. 9 
4 
| ‘ 
= 


q 


M5, 
hness 


ug 


Ro 


PIPE MATERIAL (new) FEE L. 
brass, lead,glass, transite 5x 
comm. steel, wrought iron 15x 


cast iron (avg) © 
riveted steel to 


2 4681, 


the modified Moody di: diagram obtain the m and the fric-_ 


a . Fora velocity c consistent with the Reynolds number used in 1 Step 1« of zo 


culate the ert in the pipe friction factor 


£ 


10 for Ko = 


q 
a 
a 
wea 
4 
> 


This gives values for K, and 1 m in the equation he = (= KoQ™. Because Ky and — 
m “may change with Reynolds number and hence with the equa- 


as for conditions approaching -pipe flow. Examples in the 


Appendix show in detail the method of obtaining exponential fi formulas from 
the modified ‘diagram. 


eq pipe concept is commonly used to analyze the flow for pipes” 
in series or in parallel. This may be done by replacing the individual pipes ¥ 


_ diameter but whose lengths : are chosen to give loss characteristics equivalent 
to the original pipes. Two pipes are thus considered equivalent if they hav 


equal head loss at equal rates of flow. ‘This: be expressed in terms 
A the ex nential equations as follows: oe 
eres * 


. For pipe } No. 2 2 let hy = 


In the normal use the the pipe concept the ex the exponents m4 and are 


assumed equal when a particular exponential equation is selected for the 


problem. Then if pipe No. 3 is to be made equivalent to pipe No. a the length 
of pipe No. 2 is so chosen that K.. is equal to 


_ The previous discussion of exponential formulas makes it clear that the — a 


ex _ exponents m may not be equal for all pipes ina system and if this is true 
2 there is no one pipe which can be made equivalent to the various pipes in the 


system. Equivalence may be achieved fox for one rate of discharge, but not for 
a broad range of discharge. 


_ The effect of a difference in the exponent m may be illustrated by an ex-_ 


~ ample in which a 12- inch diameter steel pipe is compared with a 6-inch . _ 
7 i diameter cast iron pipe. A 100-foot length of 6- inch cast iron pipe having a 
4 resistance relation hg = = - 1. 88ql. ‘95 may | be cummenat with 3660 ft. of 12-inch © 
steel pipe having a resistance relation hg = 1. 99Q1. 88. The length of 12- “inch 
4 pipe was chosen so that the two pipes would be “equivalent,” that is the head Ss 
loss for the two pipes would be equal, at a discharge of 2.25 cfs. The — 
oy a this comparison | are shown in Table II. It is apparent that while the pipes — 
in 
a the steel pipe is 3 12 per cent greater than that in n the 6 -inch 1 cast ir iron. n pipe, aa 
- while at a flow of 10 cfs the head loss in the 12-inch steel pipe is 11 per cent _ 
less than that in ter 6p -inch cast iron pipe. This demonstrates that when the 


conditions falling in the transition range between smooth-pipe and rough-pipe 
| Condit ng in the tran ge pip gh-pipe 
on the general resistance diagram, the results from the exponential | 
= - equation agree with those obtained from the resistance diagram within 5 per _ ae 
=a cent for a twenty-fold variation in Q. The agreement will be closer for con- Siu 
— 
> 
— 
— 
4 for Of Calculation with several nines aii of ti 
= 
a 
| 
q 
Although there are some Situations where the variation in the exponent 
is significant, in many applications of resistance formulas for design purposes 


e ect e e e a 
of 100ft. of 6-inch Ratio 
12-inch Steel Pi he, Cast Iron Fire iy 


0.49 


= 
- includes pipes with a broad range of relative roughness and a broad range of 


‘Reynolds numbers, the conventional “equivalent pipe” concept may not be © 
¥ adequate to explain the operation of the system. However, its performence | 


aaee could be understood in terms of the applicable resistance equations as deter- 

_ mined from the modified Moody diagram. In many applications of resistance | 
formulas for design purposes, the inaccuracy due to the variation in the ex- 
-- ponent m will be of little practical significance because the precision of 
-_ basic design conditions does not justify great refinement of the calculations. ‘om 


Also, in many applications it will be found that the m will | be practi - 


the same for all Pipes in the system. q 


xponential Equations ; Complex plex Pipe Systems 


In complex systems, pipes may occur in series, in parallel, or in 1 eet- 4 
works. The exponential forms of the resistance relations are especially con- 
venient for solving problems of these types. = 


— 


ae Consider n pipes in series, each with the ‘corresponding resistance inate 


— 
=a 
_ 
t 
= 
| 


= 


‘Clearly, the total ena loss for the n n | pipes in rene 


If it is possible to choose emenenibis relations s such that t the actual lia 
of each pipe can be closely approximated using the same ex 


then the t total hea head loss reduces to the simple form 


the quantity 


that the total head lo loss could b as 


~ method and subject to the same limitations referred to previously. “However, 
the resistance characteristics of each pipe is characterized by the value of © d 
Ko rather than. by some fictitious length or diameter. Equation (10) permits 
the direct solution for either the total head loss or the discharge. = 
‘if pipes: are arranged in n parallel 80 that the head loss is the same for 


> 
may be thought of as an equivalent resistance 
> 
> ig 
q 
7 
: | 
> 


in which hg, - hg. Again if it is to choose an exponential 


givi 


of parallel pipes may then be -expre ssed as 


si: 


the flow any one pipe in is desired it may be found 
equating the head loss in the particular pipe to the head loss through oe po 


pipe : systems may, for the purpose calculations, 
sited as combinations of parallel units and series units so > that the — 
i system can be analyzed by repeated application of the relationships for ai 
parallel pipes and series pipes. Thus if the Kp values are found for each 

in an exponential relationship with a constant exponent, equations (10), 
of (1), (16), and (17) will permit a complete solution in terms of the Ky values — 

ate i complex pipe system which includes some pipes which are ¢ common | 


loops can not be considered as a combination of parallel units and 


= 
Thus an equ be defined for the system 
| 
< 
| 
= 


— 
series ‘units and must noes analyzed as a network. Although the — for 
analysis of networks will not be considered here, it is pertinent to note that 


exponential resistance relations will be most convenient for use in this type 
« 


Values of the Exponent m for Water Mains 


_ Although the Hazen- -Williams equation which v uses an exponen nt m ‘of 1.85 .. 


test results appear to indicate that this is more representative of new pipe 2 


than of pipe that has been in service some time. It appears that the increased 4 


_ surface roughness on the interior of the pipe not only increases the resistance 
(giving a lower Hazen-Williams C value) but also requires a —> in — 
exponent m to a value more nearly equal to 2.0. 


Lamont(14) presents resistance data selected from more than 200 


_ of resistance measurements. In Fig. 2, the results from ninety-seven ex- I 7 
periments on new pipes of all types are plotted on the modified Moody dia- : 


smooth and rough conditions with the exponent m varying from 1.80 to 1.90. 
Fig. 3 is a similar plot for the approximately 80 records, Lamont re- 
4 ported for “old and slimy” pipes. These results clearly indicate that a larger 
value of the exponent m is required for these pipes than for new teats S For 
most of these records an exponent of nearly 2.0 is indicated. 


gram indicating that many of the results lie in the transition range between “= 


= _ Recent tests on water mains in Austin, Texas, when plotted on n the modi- 


fied Moody diagram have also indicated that values of the exponent m are as 
nearly 2.0. These results are summarized in Table III. 
‘The exponent of nearly 2.0 obtained from resistance tests on n old pipes ; 


< 
; P a ‘indicates that for pipes that have been in service for some years the resist- 


ance characteristics could better be represented by an equation of the 


7 = type or the Darcy-Weisbach equation | with a constant f than ali the 


— 


commonly used Hazen-Williams equation. 


a 


“to: a pipe for a range of Reynolds numbers. “4 
‘ =o An infinite and continuous family of exponential formulas with ‘aa exponent 
_ ranging from a minimum of about 1.70 to a maximum of 2.0 is required nd 
correctly represent all pipe flow conditions. This explains why such a multi 
oa tude of oem formulas have been presented in the literature on pipe _ 


cr 3. The commonly used exponential eihidian and the equivalent pipe concept 

plications. . With some combinations of pipe sizes, roughnesses, and flow 

‘rates, resistances may occur which cannot be explained by the common ex- x 


as based on these formulas are sufficiently accurate for many practical ap- ; 


— 


_ ponential formulas but which can be understood in terms of the concepts here . 
4 


— 
— 
3 
_ The mod ed Moody diagoram with nes O onstant exnonent mm | 
4 
) 
4 4 
4 
— 


iit 


| 


468) 


s Number = 


= 4. Exponential formulas as developed from the modified Moody digrad fora pa 
i. Be particular system of pipes may be used conveniently in the analysis of series é 
or parallel or a combination of these. 

t 


‘32 ‘5. Falganete of resistance tests on water mains as : reported herein indicate tha 
— for new pipes an exponent m from 1.80 to 1.90 is common which is in fair 

_ agreement with the Hazen - Williams formula, however, for old pipes an e 


X- 
Bi e ponent of nearly 2.0 is common which is in closer agreement to the Dare 


ie 


Coefficient in the Chezy 


a 


4 = 


r 


SS 
DSSS 


acto 


0010 


| 
vi 


— in Hazen Willams ion 
Inside diameter of pipe - ft. 

ric ion factor in Darc isbac ion. 


Gravitational constant - - - ft/sec2 


al 


ine 


= 
Coefficient in equation hg = = 


factor a coefficient i in equation hy = = K chosen to be 


— 


= 


Cement Lined 


Number 


inches 

x10° 


Cast - 
Coal Tar “18 


Cast Iron - 


Cement 0075 0.015 


Iron 


Tar Coated 008 0.013 

Cast Iron - > 


0. 


Cast Iron 
Coal Tar Lined 


Iron 


Bquivalent resistance .ce factor; 1085 for ‘series pipes; K, for 


Roughness coefficient in number of pipes | in series 
or parallel. 


Volumetric rate of fl flow - 


0.0467 5.75 0.017 0.038 2.00 
| 
2.0 
— | 
| 6 20.12 0.020 0.0008 0.0014 1.97 
| 
A 
— 
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PIPES 


Reynolds umber ad 


of piezometric head line = 


Vv Mean velocity = ft/sec 
Linear’ measure of absolute ‘roughness - ft. 
Specific weight ¢ of fiuia - -Ibs/eu ft. 


Jer 
is 


‘equation with the aid of the modified 
Find the exponential equation for 1000 fi ft. of 12- diameter 
ia ~~ pipe (e = 0. 0.00015) for velocities of the order of 4.5 ft/sec = — 


m 


the modified {= =0 0155. 


the 100 ft. of 6- inch cas cast iron pipe 


0.0018 


q 
| 
— 
3 
&§ 


ft. 


0,98 cfs” 
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CHARACTERISTICS OF RECTANGULAR THIN-PLATE rE WEIRS* 


Closure by Carl E. Kindsvater and Rolland W. Carter 


and ROLLAND W. CARTER, pa A. M. 
wl —In seniuitbia this discussion of a paper on the subject of weirs, it 
: "seems } appropriate to observe that few subjects in the literature of hydraulics 
are more revealing of the personalities of its students. 2 is indeed remarka- | 


ble that a comparatively simple physical phenomenon has inspired such an < 


4 astounding number of printed words, formulas, laboratory investigations, and 
_ disagreements. —Overa period of 250 years, nearly every name in the history 
of hydraulics is identified with contributions to the great store of weir litera- A. 
ture. . Itis not amazing, therefore, that the most common ingredient in this +f 
ma more discouraging than the repetition which characterizes sweir 
literature is | the fact that some of the best thought on the subject is buried in 
- the avalanche of printed words. . On the » other hand, the substantial area of 
ignorance regarding weirs which still exists is often dismissed with the casu- z 


__ al implication | that | the solution t to the problem is contained in ina given formula 
—which is then quoted without reference to the qualifications which were pre- 
ha Pride and prejudice seem to have played a large part in the drama of the 
simple \ weir.  Counter- -claims of accuracy and scientific truth have ben 4 
hurled across battle lines drawn up behind purely empirical formulas based a | 
on limited, often questionable experimental data. The futility of one such a a? 


=, which already occupies many pages in the Transactions, is revealed , 


_ in Fig. 4(a). Here it is clearly demonstrated that, regardless of the quality sd 


_ Of Bazin’s data, the form of f discharge equation used by Bazin and others is 
- Returning to the present, the writers are grateful for the constructive dis- 
_ cussions of their paper which were submitted by Messrs. Sarpkaya, Kolupaila, 
, Powell, Paull, Oki, and Carstens. It is regretted t that several correspondents, — 
‘ including s some respected European hydraulicians, were unable to submit 


their discussions within the time allotted for discussion. 


Proc. Paper 1453, December, 1957, by Carl E. Kindsvater and Rolland W. 


., School of Civ. 
. Regents ‘Prof. School of Civ. Eng. Georgia Inst. of Technology, Atlanta, — 
Ga. and Consultant, Geological Survey, U. S. Dept. of Interior, Washington, 


 Chf., Research Section, Surface Water Branch, Water ‘Resources 
Geological ‘Ss. of Interior, Washington, D. 
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of empirical equations. ‘They are aware the dangers O1 extrapolaung 
q empirical relations 
7 


experiments f1 from which they w were “Persons publish empirical 
4 ‘- _ equations are obliged to state very clearly the purpose and limitations of | 
_ their work. It is the reader’s obligation to observe these limitations in th * 


a 5 The writers have made it quite clear that they were concerned with the = 


verification of the formula was limited to “ water at ordinary 
u _ temperatures”. Furthermore, the range of values of each variable involved 


— 


in the experimental verification is explicit in the original paper. Without ~~ 
apologizing for the limitations which time and money placed on their work, | in] q 


the writers are, like Mr. Sarpkaya, hopeful that others will extend its bounda- 
ries. They are inclined to observe that, apart from the acknowledged limi- 14 
tations in their study of the fluid-property variables, their work covers Bes 
broader | range of geometric conditions than any other ‘single work known. a 


3 It should be quite clear from the original paper that values of k}, and kp © 4 


are dimensional quantities derived from experiments with a single liquid and 
a practical range of geometric > conditions. _ The general relationships between — : 
_ Kp, Kp, R and W are carefully detailed. One accustomed to using the caaie 
_ and Weber numbers to describe the relative influence of viscosity and surface a 
"tension, , respectively, would not expect ‘ “glucose, “honey and molasses” to 
a ~ exhibit unique characteristics as long as they behave as true fluids. er 
a Mr. Sarpkaya’ s Eq. (3) isa general statement of the slope-intercept form © 
: of the Ce curves (writers’ Fig. 10), which are believed to be descriptive of ‘3 
= the interrelationship of the principal geometric variables. His Eq. (2), on a 
_ the other hand, involves a dimensionless representation of the fluid-property 4 
ay ~ adjustment factors, kp and kp. It is not clear why he ‘presumes that an ad- 
ditional fluid- =property adjustment should be made to the relationship between — 
the geometric ratios. _ Indeed, Eq. (5) appears to hold little promise of simpli- | 


fication or clarification. And, while it is considerably n more complex than 


= 


algebraic manipulations and ‘the ‘unlikely. assumption that kp is 4 
independent of surface tension do not lend credence to Mr. Sarpkaya’ 'S pro rO- 
= tal _ Mr. Kolupaila incorrectly attributes certain objectives to pa writers (e.g. nay 
4 a search for higher precision [than + 2%] was the purpose of the new investi- 
gation”) and their sponsors (e.g. , “The sponsor. expected formula 
_ correct to + 5%”) He is also incorrect in attributing to the writers’ thoughts am 
_ which are not to be found expressed in the original paper (e.g., “Their opinion — 
7 _is based on ‘usually ass assumed values of & of about 1”, and “the conclusion of 
‘the Walchensee tests was ; expressed as being ur unfavorable by the a authors” ). 
_ The writers agree with Mr. Kolupaila that the flow pattern for rectangular — 
- : weirs might be subject to complete mathematical analysis if “a certain law of 
a streamlines” were known. ¥ Unfortunately, | the solution is hardly forthcoming a 
. ee in view of his subsequent qualification that such a law remains to be “ ser 
_ tablished or assumed” along with other “small correction» 
factors”. 


Mr. Kolupaila ‘makes reference to ‘many on weirs, 

_ most of which had been read and evaluated by the writers. Failure to refer to 
certain of these publications in the original paper is evidence that they were ae 
Ps not considered of particular value. It is not clear why Mr. Kolupaila refers — : 2 


me 
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qgictory, Dut, in general, ey appear to substantiate the writers’ conclusions 
| | 


"DISCUSSION 


to the foreign publications without identifying and supporting his « own prefer- 


ence among the diverse opinions and formulas proposed therein . A casual a 
_: of recent European reference works indicates that few of the ‘papers 
referred to have been recognized as substantial contributions in their own _ 

homelands. Finally, an examination of the paper which describes the Se 
-Walchensee tests will show that Mr. Kolupaila, not the writers, must have 
misinterpreted the conclusions drawn by Kirschmer and Esterer. The follow- 
ing free translation of their conclusions the conclusions 


“The results | of the v weir measurements substantiate the the criticism of 


reliance upon one of the rational formulae has not been n satisfactori-_ 


Mr. Powell emphasized the unreliability of coefficients derived from ex-_ 
nd periments on physically small values of P. _ The writers have acknowledged — 
en 7 q this uncertainty as well as the uncertainty of the results of tests made 
very small values of b and h in their investigation. 
we Mr. Powell’s suggested method of correcting for the e distribution of veloci- a 
ty in the approach channel is a very worthwhile and practical contribution. — 
The results which he obtained from an analysis of the Schoder and Turner a 
data are encouraging. > It is hoped that additional tests made to define the “r?_ 
_ function will refine the method. _ Thus, the weir which is unavoidably located 
in a channel which produces a non- uniform velocity distribution might still 
be » salvaged as a satisfactory meter. As pointed out by Mr. Powell, the ratio 
2 = Va/Vpb is essentially a characteristic of the installation which, like the 
_ dimensions of the weir, need be determined only once. It follows that his — 
ni of adjusting the coefficient of discharge is completely practical. It 
seems reasonable to assume that r is a function of h alone. Mr. Powell’s me 
willingness to settle for “errors of the order of one percent” isquite 
The writers | appreciate Mr. Paull’s kind cémarks. They regret. that they 
_ cannot reciprocate with an offering of dependable information regarding the 4 
_ inaccuracies to be expected under a full variety of applications. Disregarding 
such very real problems as variations in velocity distribution and ‘ — 
standard” geometric conditions, the order of magnitude of the errors to be 4 
are the scatter of the the graphs ehow- 


i proposes to demonstrate a 
lationship between C uit “red In their original paper the writers state, “The 
independent influence of this ratio (b/h) is believed to be negligible over the i | 
= practical range of the other variables. - An earlier investigation at the aaa ifs 
Georgia Institute of Technology supports this conclusion. ... A few recorded _ 
ea efforts to incorporate the b/h ratio in discharge formulas are believed to be a 
based on misinterpretations of influences related separately to the magnitudes _ 
The writers believe that this ‘statement is actually : supported by Mr. Oki’ Ss 
- Fig. 1. Indeed, closer study will reveal that the slope of the curves on the - 


figure shows the ameat of b, ‘whereas the vertical distance between the curves | 
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a between Ce (= C+ 2/3 h: for various 18 constant i 
values of b. However, if this is done, Mr. Oki’ sdata appeartoshowacon- f ) 


_ siderably larger influence due to both h and b than has been ton aapE by the © 


The writers are especially.puzzled by some of the very low values of Cg - 


shown on Mr. Oki’ s ‘Fig. 1. of Ca for than 0. 5, for ex- 


- the conclusion, based on their own work, that the h/d 3 ratio is not believed to 
_ be significant. In support of this conclusion they , submit the observation that 


the data eel in the original paper cover a range of values of h/b from 
Mr. Carstens lists. ‘several reasons for the difference in ‘the results ob- - [en 
tained by applying the writers’ method of analysis to the experimental data . 
obtained by different investigators. 7 The writers agree with his general diag - 
nosis, but they fail to see the logic in his conclusion that the maximum -experi-_ ; 
"mental error in one set of data is an explanation of the difference between _ q ‘ 
urves s which \ were drawn to represent mean v values for different sets of data. — —_ 
‘Such an unusual conclusion ignores, for example, the ‘systematic variation 
_ with velocity distribution, weir crest condition, etc., shown by Schoder and — = : 
Turner and others. To make their ove position clear, ihe writers goon — 


experiments”, a: , as the major ¢ cause e of the difference between the results of Ta 
different experiments. A large part of the tremendous store of 


Mr. Carstens’ light dismissal of the “small differences” various | 
—— ignores the purpose of most of the e work done in 1 this field since the _ 


“evidence rather than theoretical proof of the merit of the orifice analogy. 
_ The writers are pleased to endorse his conclusion that too many people write 
empirical equations (or - equations that are too o empirical 1), but they fail to see _ 
the relationship between this popular thesis and the statement quoted from 4 —_ 
_ the original paper which preceded his pronouncement. — They are also moved — 
to observe that the weir which is the subject of Mr. -Carstens’ attention (i. e., 


- b/B = 1.0) is s only a limiting case of the problem under investigation, that 
values of h/P greater than 2.0 are impractical for reasonably accurate flow 
measurement, 2 , and that the fluid-property effects associated with small va values — 
of h and b cannot be ignored for all practical pt purposes. 
In general, the discussers did not challenge the writers’ candid disrespect 
“for the traditional “theoretical” equation and the orifice analogy. However, 
most of them appeared to favor the use of a coefficient Cg which is equal to 
- Es the writers’ C divided by a quantity (2/3 y2¢ ) derived from the orifice | can 
analogy. The principal advantage attributed to Cg is that it is dimensionless. _ 
et The writers sare staunch proponents of dimensionless ratios. A review of 
their publications will subsiantiate this statement. . As engineers as well as 
q Ee) researchers, however, they strive to be e practical. In the present instance 
q _ they question the logic of straining to make the coefficient of discharge di- i 
‘mensionless when the equation of discharge contains other terms which are — 
7 dimensional. In the writers’ equation, for example, evaluation of the discharge 
: depends ¢ on | knowledge of kh and Kp, quantities which are e absolute measures of | 


— 


4 
7 
bee unique ‘in the literature. It is implied, of course, that the 
— 
a 
— 
— 
a 
P 


influences of tension. Similarly, in Eqs. 
(18) (Bazin), (19), (Rehbock), and (20) (S.1.A.), dimensional terms which are _ 
7 related to the magnitude of h restrict the application of those equations to 
It has been observed that the quantity 2/3 12g is is ‘obtained from the deri- 
a vation based on the orifice analogy. . Bazin and others have at one time pro- — J 
= equations in which 2/3 ¥ 2g was replaced by ye or Ve. - In these alter-_ 
nate forms of the discharge equation the coefficient is quite dimensionless, _ - | 

but the writers are loath to use them because of the almost exclusive division © . 

g preference between C (widely preferred by engineers in practice) and Cg a 

(preferred by some researchers) as defined herein. While they find it diffi- . 

‘ searcher finds it convenient to use C, to correlate meagre data from differ-_ ( 
ent sources. However, this reason does not justify the preservation oftheC, © 
form of the discharge > equation for practical use. % It is hoped that this brief — _ 
discussion of the relatively insignificant difference of opinion regarding 
Cq and Coe ; does not not t detract f1 from the purpose | of the siege paper. — oe 


“a Corrections. ~on page 5 of Paper 145: 


cult to get excited about the choice, the writers decided to use C mainly be- a 
- _ cause it concedes to the practicing (American) engineer the margin of ‘of con 
writers freely : admit the advantage of C than Cq. See Eq. Eq. 
c 
(15) for certain purposes. A reason seldom mentioned, for example, is that — 


Ce. is remarkably similar for different : shapes of notch weirs, Thus, the 


7 a On oa 1690-24, line 5, there should be no parenthesis after the 3.23. . 
page 1690- -25, eighth line below Eq. (F), “which Eq. should be 


“ 


2, The term V2g in the third line of the fourth ‘paragraph on page 1690- “22 
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TC TOTAL SEDIMENT LOAD OF STREAMS? 


EMMETT M. LAURSEN, LAM. ASCE. —At first reading the three ways of 
‘solving the problem of total sediment load ‘enumerated by Messrs. Garde and 
Albertson seem wise, considered and reasonable. Repeated reading, however, 
leads one to doubt, and finally reject, most of what is said or implied in this .- 
portion of their discussion. Not enough is known about fluid flow to theoreti- 
cally describe in detail the state of affairs at a complexly rough fixed bounda- 
Ty. How then does one even begin a “theoretical” evaluation of the total sedi- E 
ment transportation rate? The total load is undoubtedly the sum of the sus- 
pended and the bed load but which bed load equation does one choose and how = 
does one evaluate the reference concentration of the suspended load except by 


arbitrary assumption it (For a detailed discussion of the Einstein yocesare “4 


and the modified Einstein procedure the reader may be interested in Refer- i 


ence (6)). This leaves dimensional analysis and intuition as possible methods. i 
‘Dimensional analysis can only lead to the grouping of chosen variables. Ergo, 
"intuition must reign supreme. Unfortunately, this state of affairsis alltoo 

: true. _ Yet intuition is not entirely free, ‘but must be bound by what is known. | 


ae presented in the discussion, the relationship | for the total sediment load 


. a by intuition using dimensional analysis to form the dimensionless pa-_ 
and plotting and curve fitting to obtain the relationship between the 


— transport equations. However, its similarity to most other equations _ 
ie becomes obvious if it is rewritten. | Using the symbols of the original paper = a 
(D = yo, V, = ¥70/P) 


can be 


‘approximated by the straight line | a 


further transformed through the use of the Felation ds ct qandthe 


where the es has to absorb the various factors which 


are constant or almost so. Similarly and the critical term 


. Proc. Paper 1530, 1958, by Emmett M. Laursen, 


4 
7 
| 
a 
4 


> 


all that in some of these equations a part of the effect of the sediment size 
Mr. Bogardi’s interesting contribution, like that of Messrs. Garde and 
_ Albertson, is more a reporting of similar work than a discussion of the con-_ 
cepts ; of the. original paper. ‘The relationships presented delineating the 
configuration could be very useful if confirmed over a wider range of con- _ 
writer hopes that in Mr. Bogardi’s projected paper, the rationale 
| of at least the parameters, and preferably also the relations, will be dis- By 
cussed at greater length. The use of the total tractive force (or the shear i : 


discussion of Mr. Bondurant is s of the type that had 
~ hoped for—the application of the proposed relationships to field measure- _ 
ments. Considering the change in scale and general environment between 
4 laboratory flumes, which provided the data upon which the relationships 
“i built, and the Missouri and Arkansas rivers, the factor of four, or the shift 
inthe {(U,/w) curve, is not too disheartening. The discrepancy could be 
either’ to this gross change in scale or to deficiencies in the | proposed oo mm 
lationships, or possibly to both reasons. The fact that the original curve  de- “a 
scribed the Niobrara measurements so well and that the shift of the curve for 
the Missouri river data from Omaha was less than for the data at Kansas City — 
7 would indicate that some factor ouch as the Reynolds number might be useful 


; “the development of such a stop- gap modification and thereby enhance nell 
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DESIGN AND THE McILROY NETWORK ANALYZER 


he Closure by M. B. McPherson and Vv. Radziul > 


_ M. B. McPHERSON, 1 A. M. ASCE and J. v. RADZIUL. —The = 
mentary experiences, observations and comments contributed by Messrs. 

= and Branscome, Mr. Cole, and Mr. Lomax, 3 are sincerely appreciat- 

_ The need expressed in the paper for pittitetion ofa proven standard lard set 
of instructions for water distribution network analysis for a stored- -program > 


_ type o of computer has been | satisfied by Messrs. _ Graves and Branscome. “td 


* _ Their paper and the authors’ appeared in the Proceedings - simultaneously a 


Hamblen(1- offered a program suited to calculations for 
larger networks. - With the publication of these programs there is no longer 
a valid excuse for rough approximation by means of extensive _ shortcuts in f 
the design or investigation of complex water distribution networks by consult- — 


have not been developed,” and illustrates some of these potentials 

_by itemizing several unusual problems he has analyzed. The authors are in 7 

_ full agreement with his statement: “The combination of the analyzer and the a 
computer (at his institution) will provide complementary rather than competi- 

tive services,” but also foresee valuable collateral uses for both devices in - 


distribution system analysis by the larger utilities. — ee wii of peak day 


itatio 
raves and 


a 3 in performing the computations. Without equalizing storage, the demand - 


the peak hour of the maximum day normally constitutes the c critical > 
condition for large networks. With storage, the minimum hour of the m —_ 
_ mum day demand is apt to be more critical because the storage is i ae 
Dd being replenished at a high rate at that time and because the flow to storage " 


Bend traverses a long distance across the network. The design of ae 


‘piping improvements is a trial and error procedure. When estimated future — 
demands (always taken at rates greater than current demands) are being 
_ studied, the authors use an initial arterial network which is usually | a moder - 
ate modification of existing piping but which will probably yield e: excessive 


~ head losses. They have been able to reach a more satisfactory final vei 


- Proc. Paper 1588, April, 1958, by M. B. McPherson and J. Radziul. 
ics (1. Research Engr., Philadelphia Dept., Philadelphia, Pa. 
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arrangement more in with 2 a network 
_ which would undoubtedly prove to be more than adequate. Hence the statement — 
in the paper: needed system 1 revisions are evident at inflow rates well 
below the design rates; obviously needed piping improvements are then made - 
prior to any formal or complete test runs.” 
_ Messrs. Graves and Branscome correctly note, with reference | to Fig. 6 - 
of the paper, that the statement: “A revised storage rate requires a new ie 
storage curve,” is not clear. ‘The balancing of demands, storage rates and - sed 
_ near-uniform pumping rates is difficult to describe without : illustration. How- 
ever, a detailed example of balancing procedures using hour- -by- hour el 
- for the City of Philadelphia’s Belmont High Service District, -_ a — 
station and one storage site, be published soon. (2) 
r stepped pumping. are extremely complicated and no 
a versal rule can be applied even between service districts of the same water 
S works. The principal benefit sought in the installation of equalizing storage oi 
ist the provision of more stable pressures and less ‘Service interruption t than - 
can be achieved by direct pumping. Unfortunately, the storage volumes need-_ 
to secure near-uniform pumping, augmentation of pumping station output 
during fires, and | outage relief, are neither a summation of the volumes need- 


_ Consider first, peak day storage. As an example, if the maximum ‘aay - 
is about of the annual average, 10% of the average day 


mum day demand is s needed for. near- -uniform pumping on 1 the 1 maximum day, 
the storage volume required for the maximum day is thus about 1.5 times the 
required for the average day. ‘This theoretical ratio ranges between 
1.3 and 2.1 for the five districts cited in Fig. 4. The e ratio is undoubtedly _ - 
_ higher in many systems across the country. Maximum day storage raises es 
~ and stabilizes critical low pressures under non-emergency operating con- 
‘ditions. ‘Benefits received on days of lesser demand might be regarded as | Si. 
‘bonuses not chargeable to the cost of the maximum day storage. If the eal 
7 : provided is not sufficient to meet peak day variations, the resultant higher — 
_ pumping rates on peak days will generally require a more expensive distri- 


- bution network than would be satisfactory for near-uniform pumping with cote } 


. what period during the > design maximum day should fire requirements 
met?_ Among other things, the N.B.F.U. rates a district on ability to 
liver the : required fire flow during a district demand equal to the maximum 
day rate. ‘Philadelphia’ s design policy is to provide a working elevated | 
a storage volume (upper 25-feet) adequate to meet future maximum day demands 
rr under hypothetically perfect uniform pumping (exceeds somewhat the quantity _ 
needed for near- -uniform pumping). An additional, Separate i increment Of storage 
for fire flow augmentation is not provided. The fire flow increment is to be 
delivered from the pumping station. Special design analyses are made to 
insure adequate fire pressures under these conditions. In large districts, 
‘ “a flow and pressure requirements for fires seldom dictate additions over and 
< above the piping needed to satisfactorily meet peak day demands. | 4 
" For what duration and for what hour of the day should outage relief be ‘pro- 7 
vided? In Philadelphia it is not ‘considered economically feasible to earmark 


a increment, or piping, for outage relief alone. 
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DISC USSION 


om the storage o ona 1 peak day will be almost inconsequential. 
», the minimum amount of storage worthy of serious consideration should be the 


e - mum, near-uniform | pumping may often be feasible for m more than half of the 
year’s demands, as inferred in Fig. 3. The annual cost of pumping will be 
greatly influenced by the large share of demands thus met by near-uniform 
pumping with a single pump. The cost of an increment of storage above the 
_ minimum to meet peak day fluctuations can be readily compared against the 
cost of additional piping which would be required to provide similar service 
“a ‘pressures by means of direct peak day pumping without storage. The more | 
_ important comparison, from a cost standpoint, between the use of minimum — 
storage with stepped pumping on the peak day and either direct or near- = 
‘ uniform pumping would be purely speculative. For stepped pumping, it is al al- 
most impossible to realistically predict the most probable pattern of peak ~ - 
day demands and the sequence and timing most likely to be followed by the 
station operator in meeting these variations, 


‘The dollar value of the more ‘reliable overall service ce afforded by equaliz- — 


Discussion by M. B. McPherson and J. V. Radziul, ASCE Proc. Paper 


13-15, September, 1958. 
Cc. by J. Hamblen, and authors’ closure, Proc. Paper 
1855, pp. 11-17, November, 1958. eta. 
McPherson, B. and Radziul “Distribution Design Consider- 


ve Jour. A \.W.W.A., pending, January or 1959, 


| However, the extra benerit OF added out Her fas Deen a deciding factor 
4 adopting some otherwise marginally justified, needed improvements. 
le 
‘ses 
y 
| 7 
eel - | ly, the cost of an increment of storage especially set aside to more adequately - — 
— meet fire demands cannot be equated in direct terms to a specific expected : 
reduction in fire insurance policy rates, 
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‘Discussions by Ven Te and J. L. H. Paulhus vin 


VEN TE CHOV CHOW, A. A. M. ASCE.- —The author has presented an interesting. 
. 


rocedure for synthesizing flood frequencies through the use of a 
formula and unit hydrograph. In essence, the lengthy derivation of the pro- pro- 
cedure may be boiled down to the effect that the flood discharge can be com- 
puted by the conventional rational formula and then reduced in the magnitude 7 : 
in crder to convert the value for the rainfall frequency to that for the runoff — 
of the same frequency. When the conventional rational formula is used, it is - 
logical to accept that the frequency of the computed d discharge is compatible — 

? with the frequency of the rainfall intensity which is used in the formula. — Ow- 
* to the detention and retention characteristics ofa watershed, the actual 


above may be by a modified rational 
"where Q my the peak discharge in cfs., K is a reduction factor, C is a runoff a 
coefficient, I is the rainfall intensity in in. per hr., and A is the drainage area 
a ‘The use of the unit- -hydrograph theory in in the author’s derivation indicates es 
that the rational formula in the form of Eq. (10) is a special case of the unit- 
hydrograph method. A general mathematical ‘proof of this point has been 
Eq. (14) is in fact given in the form of Eq. q. (13) in the. paper. ‘The author pd 
— that the value of K in both formulas may be considered as a constant, “i ie 6 
_ that i is equal to 0.44 for overland flow on turf and equal to 0.77 Sorastural 2 
areas, sewered areas, and overland flow on pavements. In other words, the — 
-_ Felationship between the rainfall and runoff for a given frequency is assumed 
constant. It is believed that the author has accepted that rains of the same 
duration and frequency occur with associated runoff producing conditions that re | 
vary from time to time on antecedent factors. It may be then 


ao: truthfulness of this reasoning can only be verified by the use of actual data. 
: Although the author has shown several actual examples for this pur se, ‘it. 


Prot of Ryd 1808, October, 1958, by Franklin F F. — er 

1 Hydr. ‘Eng., , Univ. of Dlinois , Urbana, nl. 
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would be interesting to s see > how this reasoning would apply 


_ areas which are different from those investigated by the author. — ee 
generally recognized that two major factors that lead to ‘runoff 
4s the storm rainfall, a climatic factor, and the antecedent conditions, a physio- 


oe graphic | factor. The probability occurrences of rainfall and runoff are well- 


a known, while the probability occurrence of the physiographic factor is sa 


_ subject that has not yet been explored. From a probability viewpoint, how- 
ever, the physiographic factor is a statistical variable just as the rainfall and 
runoff are statistic variables. (2, 3) The frequency of ‘runoff is the product of 


Wi. ther this % 
‘is s true or not : remains to be seen from future verification 1 by d data of various 


ff 
- * _ James Edward Nash, Determining Runoff from Rainfall, Proc. Instn. Civ. a 
_Engrs., Vol. 10, pp. 163-184, June 1958. j 


a “a ah Ven Te Chow, Hydrologic Studies of Floods in the United States, Publi- 


cation 4 Symposia, Darcy, Inter. Assn, Hydrol.,  TUGG, | PP. 170, 1956, 


“the flood- frequency problem has been developed by Mr. Snyder. ‘His a 
cedure appears suitable for deriving flood-frequency data for very small 
areas, especially if they have a large percentage | of impervious s surface, , PFO 
vided representative rainfall-frequency data are available. There are, how- 
ever, several weak points in its development that make Snyder’s approach — 

less adaptable to larger | basins. In general, when representative 

a frequency data are used, the weaknesses would indicate a tendency for the 
procedure to yield results somewhat higher than might reasonably be ex- 


\ ae Perhaps the biggest drawback to the successful application of Snyder’ s pro- 
- cedure, even when properly developed, is the lack of means for selecting es 
7 - most representative rainfall - -frequency data to use. Snyder advises that vari- 
. ation in annual precipitation ‘must be considered in the selection or adjustment 
of the point rainfall record to be used but omits details. This practice does 
_ not necessarily produce reliable results. Table 1 shows the different flood- — 7 
_ frequency values that would | be obtained for | Wills Creek near Cumberland, — 
Md., by applying Snyder’, s procedure to rainfall - -frequency data, unadjusted a: 
except for partial- -duration series, for the five nearest surrounding stations _ A F 
for which that type of data are already available. Pittsburgh and Elkins are 
- each about 75 miles from Cumberland and the other three stations are between 
5 miles away. The two nearest stations yield widely different sets — 
“of results, with flood- “frequency values based on Elkins reinfall being at about 


Staff Hydrologist, U. Weather Bureau, D. 
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DISCUSSION. 


FLooD- FREQUENCY DATA. (100CFS) BASED 0 oN 5 STATIONS 


50 per cent higher those ‘Pittsburgh. is repre- 
P sentative for the basin? Noting also that Baltimore yields results over 25 
per cent higher than those of Harrisburg, with normal annual precipitation of = 
42,59 and 36.01 in. Tespectively, which of those two stations, or what kind of 


‘ adjustment, would one use in developing a synthetic flood- frequency. curve for 
7 _ The results of Table 1 would be Rinened very little by adjustment of basic — 
-frequency data for differences in mean annual precipitation. Stations 
7" _ in and near Wills Creek Basin indicate a mean annual precipitation between > 
a : _ Pittsburgh’s and Elkins’. The differences would not provide a suitable basis mn 
for an adjustment yielding flood-frequency \ values approximating the observed. 
ae = is possible that neither Elkins nor Pittsburgh yields results comparable to 
observed because Snyder’s rainfall-runoff relation is not representative 


a for Wills Creek. There isa very small chance that it c could be considering — 


— of the more serious shortcomings of the procedure is the assumption 
_ that the rainfall intensity for a given return period will produce a peak-— 
_ discharge value having the same return period. | . For example, the 10-year 
rainfall intensity is assumed to produce the 10-year flood. In this respect, 
Snyder’: s procedure is no better than formulas such as the ed ‘rational, 


es Annual maximum peak discharges often occur in connection tin rainfall 
:- intensities less than the annual maximum i intensities. The reason is that : 


many regions the highest intensities occur in the seasons when soil conditions _ 4 
4 are not generally favorable to high runoff rates. For example, Snyder would 
presumably use rainfall intensities for Washington, D. C., to develop | flood - SS 
frequency data for the Rappahannock River near Fredericksburg, Va., which © 
: has a drainage a area of 1599 square miles anda computed concentration time 
r (T,) of 24.8 hours if it is assumed that the impervious area of the basin _ 
amounts to less than 1 per ‘cent, The rainfall-frequency data he would 
“= were obtained by the Gumbel method of an analysis utilizing the annual maximum 
values for the 57-year period 1896-97, 1899-1953. The results were then ad- . 
7 J justed by Snyder for a partial duration series and converted from a - 
to amounts to yield his Fig. 2, - Appendix 1. ‘Comparison c of the ‘monthly distri 
= _ bution of the annual floods of the Rappahannock for the total 45-year period > 
(1913- 57) of peak discharge record with that of the annual maximum 24- hour — Ue 
rainfall intensities in the Gumbel analysis shows a relatively poor 


Zz It is readily seen thet: the 3- month period yielding the greatest percentage 


of annual fl floods: is March- ~May with 37. 9 per cent. Only 15.8 per cent of the Fi. 

: Berg maximum 24-hour rainfalls occurred in the same period. On the other = 
hand, while the 3-month period, July-September, yielded 57.8 per cent of the — a 


annual maximum 24-hour only 26. 6 per cent of the annual 
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2. MONTHLY DIS DISTRIBUTION OF ANNUAL MAXIMUM FLOODS: AND 


Ann. Floods nn. Max. 24-hr 


i the above for a ‘small basin for which point rainfall ae be more repre- 

sentative, it appeared advisable to make a test. Unfortunately, there are no 
"streamflow records of adequate length for any small basin in the area of ne 

interest. - Consequently, the writer synthesized a record of annual floods for __ 

the Little Falls Branch near Bethesda, Md. lane area 4.1 square miles, oa 

using: the procedure > outlined in Proc. Paper 1451.\4/ Hourly rainfall data al a _ 

Baltimore, Md., were used in synthesizing the annual flood record as they _ ' 


‘ Bearing in mind that the writer was, at this weted. concerned only with the 
comparison of the” monthly distribution of annual floods and annual maximum a 
rainfall values, it should be noted that the accuracy of the 1-hour unit hydro- 

is graph used in synthesizing the annual flood record had little bearing on their 
“8 monthly distribution. The reason : for this is that the annual maximum 1- all 
computed runoff was almost always responsible for the annual peak discharge 
the basin, which, according to the author, has acomputed T, of 1.65 
hours, Consequently, the comparison could have been made directly by using — 


. = maximum 2-hour rainfalls at Baltimore, Md., for the 49-year period, - 
1903-51, used in the Gumbel analysis for U. S. Weather Bureau Tech. Paper — 
* 25, (2) was then determined, it was noted that annual 1 maximum 1-hour rain- 
falls had almost e exactly the same monthly | distribution the 2- hour. 
comparison of monthly distributions of annual maximum floods and 2- 
Se is given in Table 3. As was expected, the latter comparison shows 
better agreement than the first. The 3-month period, July-September, has 
= the greatest percentage of annual floods and annual maximum 2- nour rainfalls 
= 
spectively. Actually, of the 49 values of observed annual | maximum 2- -hour 
wed rainfalls, 30, or 61 per cent, produced annual maximum synthesized peak dis- eo 
of “charges; other 19 synthesized annual maximum peak discharges resulted a 


from lesser rainfalls occurring with more favorable flood-producing soil 


| 
: 


e 3. MONTHLY (OF ANNUAL MAXIMUM FLOODS AND 2- HOUR 


should tend to yield flood- values th that are too high. 


fortunately, only one of the basins the author used has a streamflow record ~ 


_ of sufficient length to yield a reliable test. The Rappahannock River near . 
Ica gy Va., has a 45-year record (1913-57) of annual peak dis- 


Table FLOOD-FRE ~FREQUENCY RAPPAHANNOCK NEAR FREDERICKSBURG, | VA, 
n method j 


Inv view the author’ s a ~synthesizing procedure 


4 such as described in Proc. Paper 1451(1) would yield more uncertain results . 
than his method, comparisons were made for two other basins used by | him — 
and for which annual floods had already been synthesized. These two were © 
Goose Creek near Leesburg, , Va., drainage area 338 sq. mi., and Little Falls 
Branch near Bethesda, Md. "drainage area 4.1 sq. mi. The observed and 
synthesized annual floods w were analyzed by the Gumbel method and then a 


d- 
justed for a partial duration series. The results are listed in Table 5 for a 


comparison with the values obtained by Snyder’s method: 


by 
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FALLS. ‘BRANCH NEAR BETHESDA, MD (10 
s method (Wash. data) 167 
Snyder's method (Baltimore data) 90 ‘122 
Synthesized floods, 61 yrs. 158 
The comparative tz tabulations of flood-frequency data in Tables 4 ont 5 _ 
; cate that Snyder’s procedure yields results in fairly good agreement with the ~ 
frequency data obtained from the observed streamflow record for the ar 


4 50- and 100-year floods, whereas his values for other recurrence 
intervals may be lower or higher, suggests a . possible bias in his method. 
_ Incidentally, the writer was unable to determine how the author obtained 

__ the Gumbel values he listed for Goose Creek and Little Falls Branch. The 2 =. 

= values listed above in Table 5 can be ‘expected to be significantly higher than — 

} _ the author’s for recurrence intervals of 10 years and under because they have 
~ been adjusted for a partial duration series. However, the differences are 
i such that they cannot be explained by that alone. In both cases, the author’ Ss 
- computations appear to be at fault. In the case of Goose Creek, two annual _ 

A _ peak discharges of 45 1000 « cfs and one of 32, 800 cfs during the 26-year period — 
of record indicate that the author’s values are much too low. In the case ae ‘4 
Little Falls Branch, the writer used an additional year of record but that 2-57 
value was the minimum annual peak discharge of record. basen, 

2 there should be no appreciable difference between the frequency values from _ 
o oe 12-year record and those from the 13-year. However, the values present- 
ed by the author in his Table 2B are again appreciably lower than the writer’s: 
“a Table 5). Possible explanations are that the author obtained his Gumbel ae 
q by fitting a line visually instead of by computation and/or that he did not con- 
sider published and unpublished revisions of published annual peak ithenne. 
- Perhaps the weakest link in the development of Snyder’s procedure is the — Pe, 
ao runoff relation, which is based on monthly precipitation and runoff — ey 
data. He assumes that such a relation tends to represent storm situations and 
advances questionable reasons. _ Nevertheless, he drew his 


= = ~How much higher 2 One is inclined to conclude that the curve had to sat 
be e drawn to yield some pre- -determined 
use of a single curve to represent even average storm ‘rainfall - runoff 
cr is not sound. For example, his relation (Fig. 3, App. 1) indicates a 
© about 48 per cent runoff for 5 inches of rainfall. Presumably, it makes no 


course, contrary ‘to fact and therefore unacceptable. A duration, or 
of concentration, parameter must be included in the relation. 

— order to determine how well Snyder’ s relation fitted nena basins, a 
re relations based on monthly | precipitation- -runoff data were derived for five of ra 
_ the basins listed by him. Fig. A shows the relations for the smallest and 

largest natural basins, ‘namely, Little Falls Branch and Rappahannock River, iad 


a — Snyder’ s relation superimposed. . In both cases, as in the other cases not 
: _ shown, a seasonal effect is noted. In the case of Little Falls Branch, Snyder’s ~ g 
tion appears to to yield considerably higher than a mean curve all the 
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Figure A. MONTHLY RAINFALL-RUNOFF RELATIONS 
er ais indicate. Fig. A suggests that Snyder’s mean relation does not ba 
_ properly represent mean monthly precipitation-runoff relations, let alone — 
storm rainfall-runoff conditions, for some individual basins in the area of. 
interest. In fact, the indications are that there is a very small chance that -_ 
MMs relation would approximate average storm-runoff conditions on any basin. | 
- Another weak step in the development of Snyder’s procedure is the point- 
«area 4 relation. He assumed a random distribution of storm centers over a 
4 _ circular basin and obtained two sets of reduction factors, one for summer _ 
storms and one for winter storms. The former was assumed to apply to 
- concentration times of 4 hours and less andthe latter to 6 hours and more. 7 
_ Since frequency is involved, it would have been much better to have used the 
~ point- area relation presented in U. 3 Weather Bureau Technical Paper 29 = : 
_ to which the author refers. That relation is based on a large amount of _— ~~ 
= data of comparable frequencies from unusually dense eee 
: networks and is therefore most applicable to Snyder’ problem. nh 
a As implied above, considerable subjectivity may be involved in the a 4 
7 = of the proper point rainfall- ~frequency data for the computation of flood © 
frequencies for a particular basin. Subjectivity may also be involved in the < 
eatin of a representative value of the friction factor, n, ; and caine 
of area t drained by and the 
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percentage of natural drainage channels eliminated. ae A difference of .01 in the 
value of n may sometimes lead to about a 10- -per cent difference in the com-— 
a) puted flood-frequency values. Errors of equal magnitude may result from Brat 
incorrect evaluations of other factors used in determining Tp. 
In summary, the chief criticisms of Snyder’s procedure apply to its de-- 
velopment and not to the procedure itself. Some subjectivity was apparently | : 
involved in the development. This applies especially to the derivation of ell 
- Yrainfall- runoff relation. _ Furthermore, that relation and the point-area re- 
lation are oversimplified. Application of the procedure would also require 
some subjectivity in the selection of rainfall- data, friction factors, 
a> ; etc., to be used. In view of these shortcomings, the writer cannot agree with 
the author that his procedure would yield flood-frequency data of greater re- 


_ liability than could be obtained by synthesizing flood records. (1) However, — a 


a Snyder’s approach appears to be basically sound, and can be modified to yield 4 | 
present fo much than come | be expected from it in its 


RE FERENCES » 


1. Flood Frequencies Derived from Rainfall Data,” L. i. -Paulhus and 


J. F. Miller, Proc. Paper 1451, Jour. Hydraulics Div., ASCE, Vol. 83, Beat il 

‘Rainfall Intensity - Duration- Frequency Curves for Selected Stations in the : 
United States, s, Alaska, Hawaiian Islands and Puerto Rico,” U. S. Weather — J 

Bureau Technical Paper No. 25, Dec. 1955. 

3. “Rainfall Intensity- Frequency Regime, Part I: the Ohio 


8. 
_ Weather Bureau Technical Paper No. 29, June 1957. ris 5 ame 


i 


— 
a” 
3 
«4 
— 
fee 
I 
— 


si 


| Discussion by Ven Te Chow 


| 
VEN TE E CHow,1 A. _M. ASCE.- —The author has an extremely 
4 stimulating paper on the history, philosophy, and the goneest | situation of 4 
J hydrologic data collection and treatment at the present time. presen-_ 
tation is timely and valuable because hydrology as either a science or a tech- i> 
- nology has come of age. Owing to a general negligence of the problem, the .- 
need for an adequate basic-data program in hydrology should be > greatly 
and must be made known to all people concerned. 
_‘The use of diving rods in | the discussion of the problem in question is" 
_ interesting and appropriate. In fact, the philosophy of divining rods is not * - 
unique in the development of hydrologic while its essence is also 
true in many other scientific problems. The acceptance of witching philosophy 
the majority of the | ‘society constitutes indeed a stumbling-block to the 
progress of technology and to the development of any scientific | program. 5 
2 _ According to a recent survey made in the study at Harvard University, 
‘great of approximately 25 ore are known in existence, 


- must choose well sites under conditions where they have few external nar 

as to which of several possible sites will minimize their chances of paying — 

heavily for the water they need.” A modern basic-data program for collect- — 

ing and interpreting hydrologic information should supply sufficient “external 
aa _ guides” to those who are in need and therefore will end any attempt of _ - 


_ The gap between the data collection and the data processing and interpreta-_ 
ting is apparent from the author’s statement: “About 20,000 pages of ll 
_ data are published annually. Only a few hundred pages are devoted to inter- 
_ pretive and research results.” It is economically unjustifiable that much ef- 7 
fort and hence money have been spent on data | collection while > relatively little - 


C—O? results are produced. "¢ Although n many agencies are engaging in 


or a unified — data agency, if suowed te to be feasible, >, should do is to 


. Proc. ‘Paper 1809, October, 1958, by W. B. Langbein. 


‘ Prof. of Hydr. Eng., Univ. of Illinois, Urbana, 
; 2. _ ‘Ray Hyman and Evon Z. Vogt, 1958. Some Facts and Theories on Water- 
in States tes: GeoTimes I II (9), March 1958 1958, PP. and 


DIVINING RODS VERSUS HYDROLOGIC DATA AND RESEARCH? 
| 
all 
| 
n 
7 
wa 
is 
concerted effort or unified attempt to make the information readily available 
ae 4 aa to a researcher in the field. Many times the researcher has to spend much | ~ 
Ww energy in collecting the information from the information collected by differ - 
| 


6 March 1959 
classify and collect the information from the individual -eetbectnetn made b 


s all public and private agencies. Such classification and | collection can be os ‘ 
~ easily recorded on IBM cards and an index to the information should be issued . 
_ from time to time. The U. S. Geological Survey is currently publishing index- 

es of surface-water records in the United States. . Such i indexes should be ex- “4 
panded and extended to cover all kinds of hydrologic data by a sort of central- 
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WAVE FORCES ON SUBMERGED STRUCTURES* 


by Turgut Sarpkaya 


eros SARPKAYA, la. A. M. ASCE.— —The timeliness of this paper paper lies in. 
the fact that while | it provides design data much needed for the application i 


7 structures, it shows the method of handling the elusive problem of determin-_ 
| ing the drag in unsteady flow by means of an nw, —- careatiotic 


a This vraer s comments and criticisms will be « on testing facilities, , test st 


a wave frequencies, standing transverse waves were developed which 
produced a variation in ae height f from one side of the model to the other — 

| as much as 40 per cent. 
“height at the center of the model provided a value which differed as much as 
- five per cent from the average effective wave height. - This phenomenon was 


caused by the establishment of a resonant motion in the channel and tank, and — 
no method was found to eliminate it.” Basically the problem is that of the _ 


gravity waves in open channels. Suquet_ 
and Wallet(1) and Milne- Thomson( 2) have reported similar instabilities. 


~ Gawn(3) has demonstrated experimentally that when the frequency of the a 
oscillations of the cylindric plunger exceeds a certain well defined critical — ) 

- value, transverse waves are also set up resulting in a three dimensional _ 
gravitational oscillation. In addition to the resonant motion set up inthe 

_ model basin as described by the authors, the presence of vertical walls of the 
oma and most important of all the particular wave generating mecheninm 
se believed to be responsible for the transversal instability of waves. This 

» belief is supported by the fact that a series of plates attached perpendicularly 


to the plunger at 1 regular intervals have prevented or eliminated the ‘trans- 


bal If the authors have omalenad such anti-standing wave plates on the plunger, 
wave filters immediately at the downstream side of of the plunger, and artificial 
wave absorbers at the end of the test channel, they | might have prevented ~. 
ie. 2. This } writer believes that the authors are well justified in using Airy’ s ’ 
: _ equation for oscillatory gravity ‘waves, rather than those of Stokes or Struik, 
for the reason that the results reported(4, 5) areinbetter 


and Paper 1833, 1958, by Ernest F. Brater, McNown, 

t., Univ. of Nebraska, ibiasiion N 


q 
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with Airy’ as far as the of waves are con- 


- cerned, However, it should be noted that the wave profile is described a 
by Stokes’ or Struik’s equations. The deviation ofa curve 

the actual profile might slightly effect the AH value. 

ag 3. Computation of oscillatory wave forces or the sum of thie — ida the 

‘ inertial force on the submerged structures under consideration is based a 

a< ‘the assumption that the flow pattern is not disturbed by the presence of the a : 

am oe the flat plates the values of Cg and Cm were selected to give the 7 

correspondence e with the measured values of phase angles and the magni-_ 
= of the maximum forces. For - barge ‘models the values of phase angles Ss 
and the magnitudes of the maximum forces. For barge models the values of 
‘were computed using the measured values of inertial forces. Hence, for 
za barges there | is no way of knowing how approximate the values of Cm are, 
ae" how much justification there is for the simplified analysis. | One can only 
= that Cm is assumed to take into account whatever discrepancies may 
arise between a more realistic analysis and a simplified one, For rectang “a 
lar barge models the authors found a Crm value which varied from 1.31 to _ 

a * 76. At first these values seem to compare favorably with the theoretical an 
value of ee obtained from a two-dimensional flow analysis ignoring tl the cef- 
fects of the time dependent wake surrounding the body. If, however, one pee 
pares a plot of the rate of change of Cp, with respect to the ratio (a/b), itis 
easy to observe that for the barge model tested the value of C,, , does not vary 
more than a few per cent if only the width of the model is increased by 50%, — 

i.e., from 10 in., to 15 inches. Hence, a wake that would form m only behind the 
plate would not have any appreciable effect on the value of Cm. Since, how- fo 
ever, the sharp edges lead to separation for \ very 1 motions, s,re- 

“< sulting vortices grow into a wake and extend ess disturbed flow region also 

laterally as well as transversely. Therefore, in comparing the experimental 
— Cm values with the ones obtained from a two dimensional flow analysis, one oar, 
should take into consideration the equivalent potential flow model which is 4 


= supposed to be made baa a surface enclosing | both the body and its time - 


mote on the value of Cm plexiglass models (1/2 geometrical ratio) of ‘the: 
right rectangular p prism and the slotted barges shown in the authors’ Fig. 4 
_ were constructed in the hydrodynamics laboratory of the Engineering Me- aa 
e. chanics department of the University of Nebraska, and the C values (added 
‘mass/displaced mass) ‘were determined by means of vibratory motion, ‘For a 
is ‘full rectangular barge the value was of 0.35, and for the slotted barge 0. ks 
 Cocreapalng Cm Values are 1. 35, and 1. 65, respectively. . In the tests at 
- Nebraska, the velocities varied from 0.2 to 0. 6 in/sec ., the frequencies from. af 
20 to 60 cycles per qe and the accelerations from 40 to 400 in/sec.sec. — 
_ The above values of C,, will be used later in computing the reduction in hori- 
In what follows, this writer would like to review the recent ideas concern-_ 
_ ing virtual mass and then point out the parameters that cause marked changes 
in both Ca and Cm. . Although, considerable progress has been made on the sill 
_ experimental techniques and although the added mass coefficient has been de- 
- = for a wide range of shapes of bodies, no new physical meaning has 7 
been attached to it, and the basic ‘definition of virtual mass has remained as: = 
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virtual mass is nem: kinetic energy which suffers from the inherent werlnese 
that it assumes a constant velocity for the confined solid body, whereas all 

outward manifestations of mass should be associated with acceleration. i 
_ Therefore, the velocity must be allowed to vary and t the fact should be recog - 
bag that the added mass can yield momentum as well as as energy. ae a 

Surprisingly enough, as early as 1888, Riecke(6) recognized the essential 

“fact | that for a sphere moving with ‘constant t velocity through an infinite a 

inviscid fluid, the individual fluid particles which are pushed aside by the 

paths in its forward motion do not return to their former positions. The __ 


paths of the individual particles are not closed curves but of the shape simi-— 
wih to that shown in Fig. 1, of this paper. Hence, besides pushing the parti- — 
_ cles aside temporarily in passing, the sphere also displaces the fluid parti- 4 
cles permanently in the direction of its motion. However, the importance 


8 cylinder’ radius. Trajectories are shown 


three particles. Broken lines show the in- 
“hdl tial and final positions of particles before 
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_ mass transport and that this mass is the added mass which unveils itself 
_ The above explanation, although adding considerably to o our understanding 

of | the added mass, does not take into account the effect of the wake- or the 


or cavity as well as with their rates of change, i.e., the instantaneous magni- 
_ tude of added mass values in transient conditions depends uponthetime 
history of motion. Suggestions have been made to” approximate 
motions by parts of a sine curve and thus to Getermine a frequency parame- — 
ter _ This procedure is obviously questionable. | It is clear from the above — 
- discussion that the variation of the wake ‘geometry \ with time causes mavked ve 
changes in both Cy and C,,. Because both coefficients depend upon the state 
| development of the wake, recently McNown, (8) McNown and ne 
lated them to each other quite successfully. 
_ It should be noted, furthermore, that in an oscillatory gravity wave | wave iti Be. 
elliptical orbits 3 are Open and there i isa continuous 1 mass transport in the di- 
near r the bottom. Therefore, the value of fC, c computed from 
data ¢ also involves the effect of the inertia of wave mass transport. a ih TR 


‘volume « of fluid which replaces the s submerged structure, the e authors « could 
have used their Eq. (11), together with Eq. (4). If this is done for the case of 
barge centered at the quarter point of a wave, one obtains 


“A, ot 


or performing the simple integration, 


cos 
sh. 
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water particle orbits at a point z, and at ta a point on the still water level. — ae q 

a a Eq. ( (2) is. compared with the authors’ Eq. (16), it will be noticed that the — 

‘Second term in the parenthesis of the latter equation is missing in the former 

equation. The difference is basically due to the linearized averaging of 
horizontal accelerations along the submerged structure. . An exact analysis =) 

a made by this writer of the difference in the pressures occurring under e ot | 

crest anda trough | using Struik’s equations(10) has shown that the right side 


authors’ Eq. (18), contains more terms in the form of series 
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Cavilty=inauced mass. Fur no lor an accelerating OF decele! 5 
7 ed ric ith antaneous shane and volume of the wak 
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absolute magnitudes. This analysis is not. given 
partly because it is irrelevant t to this discussion and | partly — it is 
5 The authors have found that the reduction in displaced mass a the 
_ barge by introducing a slot is partially offset by, “the additional disturbance — 
attributable to the odd configuration of the slotted barge . . .”, and the values — poo 
_ of Cp are greater than the ones for a full rectangular barge b by | 10 to 24 per 
for values of 
e' study, made by this writer, to determine the virtual mass of multiple — 
eae lens shaped bodies, and of circular, square, and rectangular plates | 
of various sizes amply supports only the results of the above observation, _ 
- but not. the explanation | given for its cause. Since the results of this study a 
will be reported later, only part of the results obtained for parallel square © 
- plates and for sufficiently long (large ratio of length to width) parallel rec xs 
_tangular plates will be given here. 
A ad In Figs. 2 and 3, curves with dotted lines give the ratio (added mass/ : 


y displaced mass) for full right square prism and for unslotted right rectangu-_ 


lar prism, as a function of (s/a) for constant values of (t/a). . Curves with full 
lines give the aforementioned ratio for the separated square and rectangular 
plates asa function of the same dimensionless ratios. Here, only the repre- . 
sentative curves, rather than the separate experimental shown 
for the sake of simplicity. The direction of vibratory motion | is shown on the 
An examination of the figures reveals that as the plates. are » mabedia 
. from each other but a very little distance, the added mass coefficient de- A 
i. creases slightly, i.e., the fluid between the plates behaves as if it were part 
_ of the solid full parallelepiped. _ As the plates a are set farther apart, the add- 
ed mass coefficient increases, but always remains smaller than the added ed. 
: mass coefficient of a single plate. _ Hence, when the distance between the 
_ plates is not too large in comparison with the width of the plate, there isa 
_ mutual interaction of the flow patterns around the two plates. As the plates 
_ become farther apart, ie interaction becomes negligible and each plate be- 
For the slotted barge " used by the authors t/a = 1.265, and s/a = 4 7; 
_ If the point corresponding to these parameters is located on Fig. 3, it be- ‘i 
comes immediately apparent ~ each leg of the slot behaves practically as 
an independent parallelepiped. Therefore, the reason for the partial offset of = 
the reduction in displaced mass and the resulting increase in Cm is not be- _. 
cause of the additional disturbance attributable to the odd configuration of the 
slotted barge, but on tl the contrary is because of the lack of interference of the 
two sides of the ‘slot. . For the forces in the horizontal direction, a weighted — 7m 
average for Cy, can be determined easily. For the unslotted portion oA 
Cm = 1.33, and for each leg Cm = 1.95, and the weighted average, although 
approximate for this: “shape of body, is 1.55. This value is in the neighborhood — 
of what has been found by the authors for various values of T, z/d, and H/L. 


a The fact that the length of the legs is not infinite and that the ‘slotted barge t : 


model is actually more ‘complicated is, of course, not overlooked. But, this 


mentally by this writer for the two types of barges Cus = 1.65 for the yoga a 


not alter the above reasoning. Using the Cm values determined experi- 
? barge, and Cm = 1.35 for the rectangular barge) a reduction of 7 


cent in the horizontal force is computed. Jae 
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Corrections. Eq. (), a Pi term is missing due to ar an 1 obvious 
The equations for the horizontal and vertical components 0 of the orbital © 
velocity (Eqs. (4) and (5)), and the | equations for horizontal and vertical a ac : 
celerations (Eqs. (6) and (7)), and hence, the tabulated values in Fig. 5, 
not mathematically consistent as far as the rules of differentiation are con- 
cerned and as long -as the n meaning ing attached to. theta not. changed in the 
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